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Abstract 

The behavior of a Bose-Einstein ideal gas of particles in a three dimensional 
space in the presence of a uniform field, such as gravity, and of contact inter- 
action, describing the presence of one impurity, is investigated. It is shown 
that Bose-Einstein condensation can not occur. 
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The recent experiments |lj on Bose-Einstein Condensation (BEC) in atomic vapor mag- 
netically trapped, have stimulated a new interest in the theoretical study of Bose gases [0 . 
In this Letter we propose a theoretical model, that could be relevant to simulate what hap- 
pens in the presence of some particular external fields such as a gravitational uniform field 
and a point-like (or at "zero range" , or 5-like) interaction, which describes the presence of 
impurities. We start from the classical one-particle Hamiltonian of the problem, that is 

P 2 

H = 2^ + m9X3 ' ^ 
where m and p = {pi,P2iPz) are the mass and the momentum of the particle respectively. 
The contact interaction, i.e. the presence of the impurity, will be included in the formalism 
after quantization, which can be done only under the specification of particular boundary 
conditions on the wave function satisfying the Schrodinger equation || . After the introduc- 
tion of dimensionless cylindrical coordinates, 

T* 7" Z 

X\ = — cos (ft , x 2 = — sin <f> , x 3 = — (2) 

Kj Kj Kj 

< 4> < 27r , 0<r<oo, -oo < z < oo (3) 
with the quantum gravitational length and energy defined by 

1 ( V" „ S 2 « 2 

X '^ K =(s^J ' E ' = ^T' (4) 

the eigenvalues equation e — E^> # gives us the following general form of the eigenfunc- 
tions in terms of Bessel's, Neumann's and Airy's functions 



tf, iA>e (r, z, 0) = Z^e^R^Xr)/^ , (5) 

v = z - e + A 2 , (6) 

R^Xr) = A Rjl J t (Xr) + B^NiiXr) , (7) 

Z{T i )=A z Ai(r } ) + B z Bi(r } ), (8) 

where Arj, Brj, A z , B z are normalization constants to be further determined, whilst the 
dimensionless quantum numbers are defined to be 



In order to evaluate the normalization constants we have to impose the usual orthonormality 
relations 

J±Z dzZ(z, e, X)Z(z, e>, A) = S(e - e') (10) 
/ °° dr rRi(Xr)Ri{X'r) = 6{X - A') (11) 



Concerning eq. (jTlf) , it is necessary to separate the cases with / 7^ and I = 0. This distinc- 
tion is the main point of the discussion. As a matter of fact, we have to carefully take into 
account that the particles are in a uniform field and that they feel, from a quantum mechan- 
ical point of view, another extremely localized interaction (point-like or 5-like interaction) 
usually called contact interaction describing the impurity. We know that 5-like potentials 
on the plane are not mathematically well defined in quantum mechanics. The only correct 
quantum mechanical formalism to describe them is in terms of the self-adjoint extensions of 
the Hamiltonian operator which, in the present case, is given by the differential operator of 
eq. (]l|). Actually, that differential operator is only symmetric, until we specify its domain. 
These self-adjoint extensions of the Hamiltonian are non-trivial in the case 1 = 0, when 
both the regular and the irregular wave functions at r = are allowed, in contrast with the 
case I 7^ 0, in which singularities at r = can not be accepted. It is indeed possible to see 
from the behavior of the Neumann's function for large argument @] that Br i = 0, for the 
wave-function has to be locally square integrable. Therefore we readily get 

A R)l (\) = -V\, 1^0. (12) 

In the case I = the presence of the irregular part of the wave function is allowed since it 
turns out to be locally square integrable. A straightforward calculation leads to the result 
that the normalization condition (|ll|) is true if and only if 

7rtan[7r^(ao, A)] = ln(a /A 2 ) , A > , (13) 
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where < «o is an arbitrary constant parameter. Moreover, taking eq. fllTf ) into account we 
definitely obtain 

^4j?,o(A) = —VXsmfa) , B R>0 (X) = v / Acos(7rz/) . (14) 

Now it is very important to understand the role of atQ. As we have already emphasized, we are 
dealing with a quantum-mechanical point-like interaction, which can be correctly handled 
through the formalism of the self-adjoint extensions of the Hamiltonian operator ||. The 
latter formalism entails that in order to describe the dynamics a whole one-parameter con- 
tinuous family of different quantum mechanical Hamiltonians has to be considered. The dif- 
ferent self-adjoint Hamiltonian operators are naturally labelled by the parameter ao- When 
«o = (Friedrichs' limit), the contact interaction is turned off and, from the mathematical 
point of view, it means that the irregular part of the wave-function disappears. Actually 

tan(-7rz/) = -oo =>- v = -1/2 , (15) 

i.e., the same result as in the case I 7^ 0. Let us finally evaluate the constants Az,Bz- 
Substituting eq. (Q) into eq. (|I~0D, and noting that from the asymptotic behavior of the 
Airy's function Bi(z) we have to set Bz = - in order to keep square summability on the 
positive half-line z > - we finally get 

/oo 
dzAi(z, e, X)Ai(z, e', A) = 5(e - e') . (16) 
-00 

Using the integral representation for the Airy's function || we find 

A z = l/V2^ , B z = , (17) 

the wave functions of the eigenstates with vanishing angular momentum and with a contin- 
uous degeneracy A > due to the transverse momentum becomes 

* 0iV (r, z, 0) = /^ A ^ V \ [JVb(Ar) - tan (ttz/) J (Ar)] , (18) 
Zix sec (vrz/J 

whereas the eigenf unctions of the states with non-vanishing angular momentum read 



%, x , e (r, z, 0) = J,(Ar)A%)e^ . (19) 

From the previously obtained values for the normalization constants, it follows that the 
improper eigenf unctions are normalized according to 

(^z,vI^',av> = 6w8(\ - X')5(e - e') . (20) 

Now we are ready to discuss the existence of a characteristic state which faithfully 
and uniquely specifies any given self-adjoint extension, within the above mentioned one- 
parameter continuous family of the quantum Hamiltonian operators that are allowed ac- 
cording to the general principles of quantum mechanics. To this purpose, let us start again 
from the Schrodinger stationary equation, after the replacement A — > iX. Then we obtain 

' f rf] J*(r) - AV - I 2 = , (21) 



Ri{r) dr \ dr 
1 d 2 



Z(z) -z+e+\=0. (22) 



Z(z) dz 2 ' 

The first equation is exactly the same that we had previously discussed in the case of states 
with a continuous degeneracy in the transverse momentum. Our interest relies only in the 
subspace of vanishing angular momentum because it is only in this subspace that non-trivial 
self-adjoint extensions of the Hamiltonian actually occur. For I ^ the solution of eq. fl2T|) 
is not acceptable because it does not belong to L 2 (R 3 ). Taking I = and evaluating the 
wave function of the state that characterizes the extension we get 

*o,e(r, z) = 4=^0(1 A | r)Ai(z - e - A 2 ) . (23) 

This state has to be orthogonal to any of the improper states with a continuous degeneracy 
of the transverse momentum A > - see eq. (|18]). We expect to find that the orthogonality 
takes place when the parameter A 2 of the state characterizing the self-adjoint extensions is 
exactly equal to the above introduced parameter a . Actually we require 

*o I A,e|*o, e ') = , Ve^e' , A>0. (24) 



Explicit evaluation leads to the result 

7rtan(7rz/) = ln(A 2 /A 2 ) , (25) 

which exactly corresponds to the definition of a in eq. QTBp. Consequently, orthogonality 
occurs iff A 2 = ato, and from now on we can label the state that uniquely characterizes the 
extension with the parameter «o- 

The calculation of the one-particle Heat-Kernel is the first step to get the partition 
function of the system. It is also interesting to examine the limiting cases in which the 
background fields are switched off. As a matter of fact, we have to recover the Heat-Kernel 
of the non-relativistic free particle when both of the external fields are turned off, i.e., g — ► 
and cto —>■ . To get the general form of the Heat-Kernel, we need the spectral decomposition 
of the quantum Hamiltonian operator. To this purpose, it is necessary to keep in mind 
that we have to deal both with the states with a continuous degeneracy, labelled by the 
eigenvalues of the transverse momentum, and with the single state characterizing the self- 
adjoint extension. Furthermore, it is convenient from now on to measure energies in units 
of quantum gravitational energy scale E g . As a consequence, we define the dimensionless 
Hamiltonian operators as 

E- x H(ao) = H(ao) = HK« )] + H(a ) , (26) 

where H is the part of the spectral decomposition of the Hamiltonian that involves the 
states with a continuous degeneracy of the transverse momentum, and H is the one that 
involves the characteristic state of the self-adjoint extension. The explicit form of the spectral 
decompositions can be written as 

HK«o)] = J£ de e / °° dX £, eZ |e, A, I >< e, A, l\ , (27) 
H(«o) = de e | e, a , >< e, a , | . (28) 

Let us come back to the evaluation of the diagonal kernel. For a matter of simplicity, we 
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shall calculate separately the part of the kernel corresponding to the continuous degeneracy 
in the transverse momentum, and the other one which characterizes the extension: namely, 



G(a , 7; r) = G(a , 7; r) + G(a , 7; r) , 
G(a , 7; r) = e 72 (r | exp{-^E.[u(a Q )]} | r) , 
G(a , 7; r) = e 72 (r | exp{-7H(a )} | r 



(29) 
(30) 
(31) 

where r = (r, <f), z) is the position in the dimensionless coordinate space associated to the 
system whereas 7 = f3E g and j3 = l/fc^T, with ks the Boltzmann's constant. The reason 
for the factor exp{7z} in eq.s 

(|(]HI) is that 

one has to recover translation invariance in the 
absence of contact interaction, i.e., in the limit ao — > 0. The diagonal kernel G is that one 
with a continuous degeneracy in the transverse momentum. From eq.s (|30|) and (|27D we can 
write 



/oo rco 
dee^ / d\J2\%x,e 
-00 Jo leZ 



y,z,(j))\' 



(32) 



Substituting the wave functions (PD-flU) into the expression (^) we obtain 

1 



G(a ,r,r) 



(2tt)- ./-.x 



drje^Ai 2 ^) / dAAe~ 7A 
Jo 



X 



COS ( 11V ) 



Jl{\r) - Nl{\r, 



+ sin(27ri/)Jo(Ar)AT (Ar)} 



(33) 



The integral over the Airy's function can be further elaborated and eventually set in the 
more suitable form where 



/+00 
dr, e^Ai 2 ( V ) = h + I 2 + I 3 
-00 



(34) 



with 



h 
h 



1 

18 V t I 2 



£ 3 

1 f°o JC (-7) 

,/9^ J-OO 



Kl) 



2 

3 7 



e 2 

— COS T77 , 

l€l(7 2 +C 2 ) 12 
£ 2 

Sin frr , 



^T(7 2 +? 2 ) 12 

7 



(35) 
(36) 
(37) 



E p (z,fj) being the Mittag Leffer function ||. 

Finally, taking the above integral representation into account, we can rewrite the Heat- 
Kernel G in the simpler form 

G(a , 7; r) = /(7)/(8tt 2 7) {1 - / °° dt e"* cos 2 (™) 
x J 2 (sy/t) - iV 2 (sVt)] - sin(27rz/) J (sy/i) N (sv^) } , 

s = V^g+^D ) tan (7rz/) = ln(a 7/ In t) , (38) 

where Ay is the thermal wave length. 

The evaluation of the diagonal kernel G can be done in close analogy with the previous 
case. From eq. ( |23"D with A 2 replaced by a and from eq.s ( |2H ) and (RT) we get 

G(a , 7 ;r) = Je^^(r yfa)I(>y) ■ (39) 

It is very interesting now to examine the behavior of the kernel in the limiting cases in 
which the background fields are switched off. Particularly, we look at the behavior of the 
diagonal Heat-Kernel when g — > 0, i.e., when the uniform field is off, and when oq — > 0, that 
means without point-like impurity. Finally we will examine the case when both the external 
fields are turned off and we show that there is commutativity with respect to order of the 
operations. We do emphasize that, in order to evaluate those limits, it is necessary to pass 
from the dimensionless Heat-Kernel, to the corresponding one in physical units: namely 

G{a , g, T; x) = K 3 G(a , 7; r) = 0^-) G{a , 7; r) , (40) 

where the dimensional quantities are labeled with a tilde. From the behavior for small 7, of 



the integral 1(7) ~ y (tt/t) - see eq.s ( j34j - |3~7|) - it immediately follows that 

G(g -> 0) = A^ 3 . (41) 

which shows that the above discussed specific form of the contact interaction does completely 
depend upon the presence of a non-vanishing uniform field. In order to perform the limit 
in which the point-like impurity is removed, it is important to keep in mind that ao and v 
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are strictly related by eq. fllTf), so that v goes to (—1/2) when a goes to zero. The explicit 
evaluation yields 



which exhibits manifest translation invariance as it does. From eq.s (f4T|) and (0), after 
switching off the remaining background fields we get the expected result 

lim G(a , g, T; x) = A^ 3 = lim lim G(a , g, T; x) , (43) 

g^O g^O a ->0 

i.e., the Heat-Kernel of a free non-relativistic particle. Furthermore we have shown that the 
final result is independent from the order with which the external fields are switched off. 



The partition function per unit volume of a single molecule is now to be evaluated in the 
thermodynamic limit. To this aim, let us first consider the case without contact interaction, 
i.e., «o — > 0. In this case, owing to translation invariance, we immediately get 

/2mE g \ 3/2 I X/3 E g ) 



Z (f3) = / dE p(E)e-? E 
Jo 



71 



V h 2 



(3E g 



(44) 



The corresponding density p(E) of the one-particle quantum states, i.e., the number of the 
one-particle quantum states per unit volume and within the energy interval E and E + dE, 



E > 0, can be obtained as the inverse Laplace transform of the above expression (fH). 
Explicit evaluation yields 



p(E) = ttJE, 



/2m\ 3/2 
Vh?) 



E 
E, 



-At 2 - 



E 

~En 



(45) 



As a consequence, we can write the average density of particles in the form 

(N) T 



V 



<n> 3 



°°dE^Z, 



o 



(46) 



l- Z e-0E V(l-z) 

where z = e^ M . From the above expression it immediately follows that there is no Bose- 
Einstein condensation in the presence of a uniform gravitational field, owing to the behavior 



Now, it is quite evident that we can rewrite eq. ([46]) in the suggestive form 

g x {z) , _ f» Jp p(£)-p(0)^ 



(w)r = 2^(1/3)]^ + Z io dE 1 - ze-^ 6 

+ 7(1^) ' (48) 

which shows that the lack of Bose-Einstein condensation in the presence of a uniform field 
in three spatial dimension just corresponds to the very same phenomenon for an ideal gas 
of free particles in two spatial dimensions. A little though readily drives to gather that also 
in the case «o 7^ the same conclusion holds true, namely condensation does not occur, in 
contrast with the claim of ref. Iffll As a matter of fact, if we set 



G(a , 7; r) = G(a = 0, 7; r) + AG(a , 7; r 



^7T 2 7 



AG(a ,7;r), (49) 



it can be proved f| that integration with respect to the radial dimensionless coordinates 
leads to the finite result 

| o °° dr r AG(a , 7; r) = f dt^^r . (50) 

It immediately follows that in the thermodynamic limit the only relevant term is the first 
one in the RHS of eq. (fi9"|), all the rest disappearing in that limit, i.e. no Bose-Einstein 
condensation is possible in the presence of a uniform field and of a point-like impurity. It 
should be noticed that, from the very same results of ref. M, the absence of Bose-Einstein 
condensation persists even if we describe the impurity by means of an Aharonov-Bohm 
vortex 0. 
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